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Abstract 

Supermassive black holes exist in the centers of galaxies, including Milky Way, but there is no 
compelling theory of their formation. Furthermore, observations of quasars imply that supermas- 
sive black holes have already existed at some very high redshifts, suggesting the possibility of their 
primordial origin. In a class of well-motivated models, inflationary epoch could include two or more 
periods of inflation dominated by different scalar fields. The transition between such periods of 
inflation could enhance the spectrum of density perturbations on some specific scale, which could 
lead to formation of primordial black holes with a very narrow range of masses of the order of 10^ 
solar masses. These primordial black holes could have provided the requisite seeds for the observed 
population of supermassive black holes. 



I. INTRODUCTION 



Supermassive black holes (SMBHs) with masses 10^ — IO^'^Mq reside in the centers of 
most galaxies [IH3], including Milky Way j^-Ej. Furthermore, observations of quasars (QSO) 
with redshifts z > 6 imply that SMBHs must have already existed at such high redshifts. 
Stellar explosions can produce black holes with masses up to 10 — 15AfQ [7J, but there is 
no mechanism by means of which such small objects could grow to become SMBHs [8], 
except if dark matter has a sufficient self-interaction to facilitate a rapid transfer of angular 
momentum and kinetic energy [9] (such models exist, but they are somewhat ad hoc p^ITT] ). 
The early formation of SMBHs, which is necessary to account for high-redshift quasars, 
implies that SMBHs may have preceded star formation [12]. The masses of SMBHs exhibit 
a remarkable correlation with the bulge masses of their host galaxies [21 [T2HI1] . The bulge 
mass is 10^ times larger than the black hole mass, and the proportionality holds over some 
four orders of magnitude. 

The robust proportionality between the mass of the SMBH and the mass of the bulge, and 
the early QSO activity suggest that the origin of SMBH may be primordial. A number of 
processes in the early universe could lead to formation of primordial black holes [15] . Collapse 
of primordial density fluctuations [16j usually results in a broad distribution of masses, 
unless some speciflc scale is favored, for example, by a flrst-order phase transition [17], or 
by formation of domain walls in a second-order phase transition [IB] . The idea of primordial 
origin of SMBHs was examined in Ref. [19] for the class of models that generate a broad, 
power-law spectrum of black hole masses. In such models, the required number of SMBH 
seeds can only be produced at the expense of generating too many smaller-mass black holes, 
which leads to a conflict with existing upper bounds. Furthermore, aside from observational 
bounds, the history of black hole formation and mergers favors a rather narrow mass function 
for primordial seeds. Bean and Magueijo PD] have devolved the black hole mass function 
inferred from observations and obtained the required distribution of seed masses, shown in 
Fig. [Tj This distribution, is extremely narrow, much more so than the distributions expected 
from most models of the early universe. 

However, it was recently pointed out that, in a class of hybrid double inflation models, 
a very narrow distribution of primordial black hole (PBH) masses can be produced [2TI - I27] . 
Furthermore, the masses of PBH produced by inflation can lie in a very broad range, and 
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specific models have considered masses from as small as 10~'^Mq and as large as lO^M©, 
including a particularly interesting range around O.IMq, which could be probed by gravita- 
tional lensing of massive compact objects in the Milky Way halo [231 EH]- Scalar fields with 
flat potentials are generic predictions of many theories beyond the standard model, and a 
two-stage or multi-stage inflation can actually be a generic feature of supersymmetry and 
string theory [291 13^. String theory predicts a large number of scalar fields, moduli, such as 
compactification radii, which appear as low-mass fields in the low-energy effective field theo- 
ries. Supersymmetry (in theories with or without superstring UV completion) predict a large 
number of flat directions in the potential. Some flat directions are parameterized by scalar 
fields that gain relatively small masses from supersymmetry breaking. The light scalars 
present a very difficult problem for cosmology [311 132] because their coherent oscillations 
can come to dominate the energy density of the universe very late. With the help of thermal 
inflation [331 131] , or by means of coupling the moduli fields to the inflaton [33 136] , the moduli 
problem can be ameliorated in some range of masses [37H3^ . but neither approach provides 
a complete solution. However, some scalar fields with relatively flat potentials can serve as 
inflatons and allow for a cascade of sequential stages of inflation. Each of these stages may 
be too short to account for the observed flatness and homogeneity of the universe, but they 
can be responsible for generating primordial black holes that can play the role of dark mat- 
ter [211 [39]. The moduli fields can generate density perturbations [IQ]. Given the plethora 
of hght scalar fields expected in the landscape of string vacua, a multiple-stage inflationary 
cascade appears to be a rather generic possibility. We will concentrate on the latest two 
stages of this potentially elaborate cascade. 

Since it is likely that more than one inflationary epoch could have taken place in the 
early universe, the slope of the scalar potential could have an abrupt change in transition 
from one inflation to another. The transition between inflations can enhance the spectrum 
of density perturbations on some specific scale. When these density perturbations re-enter 
the horizon, they can form a population of black holes with masses in a narrow range. 
The masses produced in these models can be as large as lO^M©, which is consistent with 
the range required for seeding SMBHs shown in Fig. [Tj This range is also consistent with 
the observational constraints. Assuming that every galaxy has a seed PBH of mass ~ 
(10'^ — 10^)Mq, the fraction of dark matter in such objects is below 10^^, which is below the 
observational upper limits [HI |12] . 
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FIG. 1. The devolved mass distribution of PBHs (on the left) that produces an acceptable evo- 
lution of mass and number density of black holes in galactic halos (right), according to Bean and 
Magueijo |20j . The data points are from Ref. |43j . The two distributions are for the limiting halo 
mass of W^Mq (dashed line) and W^^Mq (sohd hue); see discussion in Ref. [20] . 

II. PBH FORMATION IN A DOUBLE INFLATION MODEL 

Here we describe the basic idea how a double inflation model produces the density per- 
turbation spectrum with a sharp peak. After inflation an inflation field begins to oscillate. 
During this oscillating phase, if the inflation has an appropriate coupling with another field 
ip, it decays rapidly into ip through parametric resonance. Because of the band structure of 
the parametric resonance, the inflation dominantly decays into ip particles with a specific 
momentum kp. In terms of fields, the parametric resonance induces large fluctuations of 
the ip field in a very narrow wavelength range. In general, the enhanced ip fluctuations 
have a momentum of the order of the inflation mass m^, which corresponds to length scales 
smaller than the horizon size, and, therefore, they do not affect the density perturbations 
on cosmological scales. 

However, if there is a second stage of inflation following the first one, the situation changes. 
Let us suppose that the second inflation is caused by some scalar field ip. When the second 
inflation ip has a coupling with ip field, the enhanced ip fluctuations become a source of cp 
fluctuations during oscillating phase of the first inflation. As a result, ip has large fluctuations 
Since these fluctuations are exponentially extended by the second inflation, they 
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generate cosmologically relevant density perturbations whose spectrum has a sharp peak. 
This mechanism, which produces a sharp peak in the spectrum of density perturbations, 
was first considered in the supersymmetric smooth hybrid new inflation model [23] (see 
also |2ll|25l|27]) and then in Kahher moduli double inflation [H]. Evidently, the mechanism 
operates in a variety double inflation models. 

Now let us consider the relation between BH mass and the parameters of double inflation. 
The PBHs are produced when the density perturbation re-enters the horizon, k/a = H, and 
their average mass is approximately given by the horizon mass. Since the spectrum has a 
sharp peak at A; = kp, produced PBHs also have a specific mass Mbh corresponding to kp. 
Mbh and kp are related: 

Mbh - 2 X 10^^ ' , (1) 

/ n.f \ -1/2 

fcp^4.6xlO«Mpc-M . (2) 

Here kp is a comoving wavenumber. The peak wavenumber kp is of the order of the inflation 
mass, 

k 

— = am^, (3) 

'^osc 

where Oosc is the scale factor at the start of the (p oscillation after the first inflation and 
a = O(O.l). Let us define t^: and a^: as the time and the scale factor when the pivot scale 
{k^ = 0.0002 Mpc^^) leaves the horizon during the first inflation, a* is determined by 

-=Hu (4) 

where Hi is the Hubble parameter during the first inflation. Then, the number of e-folds 
A^^i* from to the end of the first inflation is 

A^'i^, ~ 14 for Mbh = IO^Mq and ln(ifi/am^) ~ 1. The pivot scale corresponds to A^^, = 
50 — 60, so the difference A^* — A^'i^, should be provided by the second inflation. 

Besides the sharp peak, the double inflation predicts some characteristic features in the 
spectrum of density perturbations. Because perturbations on large scales and those on 
small scales are produced during the first and second inflations respectively, their amplitudes 
are generally different. The amplitude on large scales is fixed by CMB observations |15], 



but the density perturbations on small scales can be much larger or smaller depending on 
models of the second inflation. In addition, the double inflation predicts a relatively large 
running of the spectral index dus/li^k where is the spectral index. In single-field slow- 
roll inflation, dng/lnk is negligible at observable scales because of smallness of slow-roll 
parameters. However, in the double inflation the observable scales exit the horizon near the 
end of the first inflation when the inflaton rolls down the potential faster, which leads to 
increase of running of the spectral index. 



III. SMOOTH HYBRID NEW INFLATION MODEL 



In this section, we discuss an example of a double inflation models which produces a 
sharp peak in the spectrum of density perturbations. The model we adopt here is the 
smooth hybrid new inflation model [23] ■ The original model was built in the framework 
of supergravity, but the essence of the inflationary dynamics is described by the following 
scalar potential reduced from the full scalar potential in supergravity: 

V = Vn + V^ + Vhn, (6) 



2 / V ^ 16M2/ 16M4 



VM = v'(l--V>')~'^vV + ^v', (8) 

v.M^,,) - (-.^ + T^)^ ? - (V + ^) .v., (9) 

where and ip are the inflation and waterfall fields of the first inflation (= smooth hybrid 
inflation [46j), ip is the inflation field of the second inflation (= new inflation |17]), n and v 
are the scales of the smooth hybrid and the new inflations {fi > v), M is some cut-off scale, 
and K, and g are coupling constants. Here we set Mp = 2.4 x 10^*^ GeV = 1. The scalar 
potential consists of three parts, Vr, and Vhn which are responsible for smooth hybrid 
inflation, new inflation and their interaction, respectively. 

For (j) > a/ jjM, the potential Vn has a local minimum at ^/^ ~ where the potential 
is described as 



This potential is dominated by the vacuum energy /i^ and smooth hybrid inflation takes 
place. The Hubble parameter of the smooth hybrid inflation is given by Hi = jj? / a/3. 
During the period of smooth hybrid inflation, the potential Vhn can be written as 

VHN = ^/iV' + /"V0yp + .... (11) 

Thus, the inflaton of the new inflation is stabilized at ~ — (t'//i)^0, which evades the initial 
value problem for new inflation. 

After inflation, the inflaton fleld and the waterfall fleld if) rolls down toward the true 
minimum = and ip = 2^/JIM and start oscillation. During the oscillating phase, fluctu- 
ations of both flelds, 5(f) and dip grow rapidly through parametric resonance due to self- and 
mutual couplings. The fastest growing mode has wave number 

-^0.3m^, (12) 



*osc 



a, 

where = ^/S/i^/M is the mass of the inflation 0. Then, the fluctuations 50 and 6i/j 
amplify the fluctuations of the inflation (p of the new inflation via the interaction represented 
by the second term in The most amplifled mode has the same wave number as that 



given by Eq. (12). After and ip decay by parametric resonance and usual perturbative 
process, the vacuum energy f ^ in Vn dominates the universe and the second stage of inflation, 
i.e. new inflation starts. As a result, the enhanced fluctuations of (p produce a sharp peak 
in the spectrum of the curvature perturbations. 

We obtain the power spectrum by numerically solving the evolution equations for linear 



perturbations of the flelds and metric,*^ The resultant power spectra of the curvature 
perturbations, V(^{k), are shown in Fig. [2j In obtaining the spectra we take /x ~ 2 x 10~^, v = 
fi/4,K, = 0.05 and g = 2 x 10^^. The solid, dashed and dashed dotted curves are for 
M = 0.52, 0.59 and 0.48, and they have narrow peaks at ~ 6 x 10^ Mpc"\ 2 x 10^ Mpc~^ 
and 2 x lO'' Mpc~^ which correspond to PBH masses ~ 5 x IO^Mq, 5 x IO^Mq and 5 x IO^Mq, 
respectively [Eq. ([T])]. The substructure of the peaks is due to resonance bands. In Fig. [2]the 
curvature perturbations on large scales [k < 1-10 Mpc~^) are produced during the smooth 
hybrid inflation and their amplitudes and spectral indices are consistent with WMAP 7-year 



In the calculation we do not include non-linear effects and back-reactions which may be important when 
the parametric resonance occurs. However, some part of back-reaction is taken into account by taking 
large decay rate for the inflation and waterfall fields of the smooth hybrid inflation. 
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FIG. 2. Power spectra of curvature perturbations for /i ~ 2 x 10~^,f = /i/4, k = 0.05 and 
g = 2 X 10^^. The solid, dashed and dashed dotted curves are for M = 0.52,0.59 and 0.48, 
and their peaks are responsible for formation of PBHs with mass ~ 5 x 10^Mq,5 x IO^Mq and 
5 X W^Mq, respectively. 



data V^{K) = (2.43 ± 0.11) x 10"^ and Us = 0.963 ± 0.014 [45J. On the other hand, the 
new inflation produce the curvature perturbations on small scales {k > 1-10 Mpc~^) and 
the amplitude is much larger. Since the observed spectrum is almost flat at A; < 1 Mpc~^, 
the new inflation should only account for the perturbations at A; > 1 Mpc^^. This sets the 
upper limit on the e-folding number of the new inflation and the mass of PBHs from Eq. Q 
which is given by Mbh ^ a few x IO^Mq. 

A more stringent constraint on the amplitude of the power spectrum comes from the CMB 
spectral distortion due to photon diffusion [50]. After the perturbations re-enter the horizon, 
fluctuations in photon-baryon fluid starts acoustic oscillation. During acoustic oscillation, 
photon diffusion dissipates fluctuations on small scales. Then the energy of the acoustic 
oscillation of the photon fluctuations is injected into the background radiation and distorts 
the CMB spectrum unless the injected energy is thermalized by photon number changing 
processes such as double Compton scattering. For a given redshift z the critical scale k^^ 
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below which the fluctuations are damped is 



k^^ ~ 2.5 X 10^(1 + z)-^/^ Mpc. (13) 

When the double Compton scattering becomes ineffective at ^dc — 2 x 10^, the critical 
scale is given by A;d(-2Dc) — 10'' Mpc~^. Thus, the photon fluctuations with k < fcdl^^oc) 
dissipate after zj)c and cause the spectral distortion which is parametrized by the chemical 
potential fi. From the COBE observation fi should be < 10~^ [51] • In our model fi is given 
by ~ j'''^ d\nkP(^{k) which amounts to O(10~^) if kp < fed- So we must require kp > k^, 
which leads to the stringent upper bound on the PBH mass as Mbh ^ 10^ Mq. Therefore, 
the present model can produce PBHs with mass 10'^~^Mq which evolve to SMBHs. 

As mentioned in the previous section, double inflation models predict a large running 
spectral index dris/dlnk. In Fig. [s] the contours for PBH mass IO^Mq and IO^Mq are 
shown in the Us-dng/dlnk plane. The spectral running is more signiflcant for larger PBH 
mass because the observable scales exit horizon at later epoch of the flrst inflation when the 
slow-roll parameters are larger. The predicted spectral running is consistent with the recent 
CMB observations by SPT [18] and ACT [49] which suggest negative spectral running. If 
the spectral index and its running are determined more precisely, it would be a good test 
for the PBH formation mechanism by double inflation. 

IV. CMB AND REIONIZATION 

Primordial black holes can absorb gas and produce ionizing X-ray radiation, which can 
have observable consequences [^ 1^ . X-rays emitted by supermassive primordial black 
holes can produce ionization and heating of gas during the "dark ages". While direct ion- 
ization is probably too weak to affect the optical depth measured by WMAP, the increase 
is signiflcant enough to enhance formation of molecular hydrogen, which is an important 
cooling agent needed for the collapse of gas clouds leading to formation of the flrst stars. 
WMAP constraints [42j allow for one supermassive primordial black hole per galaxy, as our 
model predicts. The effects of increased H2 fraction on star formation [41] can be probed 
by future 21-cm observations. 
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FIG. 3. Spectral index and the running parameter predicted in the smooth hybrid new inflation 
model for Mbh = IO^Mq and W^Mq. We also show the constraint from the combined analysis of 
WMAP and SPT data 08] 



V. CONCLUSION 

In this Letter we have considered primordial black holes (PBHs) as seeds for supermassive 
black holes (SMBHs), and we studied their formation in double inflation models. To seed 
SMBHs with masses 10^~^^Mq, PBHs should have masses as large as 10^~^Mq. We have 
investigated the smooth hybrid new inflation model as a model of double inflation and 
numerically solved the evolution of inflaton fields and their fluctuations. It was shown that 
the smooth hybrid new inflation generates a spectrum of the density perturbations with a 
sharp peak which leads to formation of PBHs with masses ~ IO^Mq. The sharp peak in the 
spectrum ensures that the resultant mass distribution for PBHs is also very narrow, which 
is consistent with the mass function inferred from observation pU] . 

We have adopted a specific model for double inflation in the present Letter. However, 
double or multiple-stage inflation is quite generic in supersymmetry and string theory, and 
inflaton decays due to parametric resonance in preheating epochs between inflationary stages 
are also common. Therefore, it can be expected that PBHs with a narrow mass distribution 
are produced in the inflationary universe. Besides the sharp peak in the density perturbation 
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spectrum, double inflation models generally predict a large running of the spectral index. 
For the smooth hybrid new inflation model we obtain —0.02 < [dng/dlnk) < —0.01, which 
can be tested in future CMB observations. 
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Initiative (WPI Initiative), MEXT, Japan. 
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